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Abstract 

We study the lubrication of fluid-immersed soft interfaces and show that elastic deformation cou- 
ples tangential and normal forces and thus generates lift. We consider materials that deform easily, 
due to either geometry {e.g a shell) or constitutive properties {e.g. a gel or a rubber), so that the 
effects of pressure and temperature on the fluid properties may be neglected. Four different system 
geometries are considered: a rigid cylinder moving parallel to a soft layer coating a rigid substrate; 
a soft cylinder moving parallel to a rigid substrate; a cylindrical shell moving parallel to a rigid 
substrate; and finally a cylindrical conforming journal bearing coated with a thin soft layer. In 
addition, for the particular case of a soft layer coating a rigid substrate we consider both elastic and 
poroelastic material responses. For all these cases we find the same generic behavior: there is an 
optimal combination of geometric and material parameters that maximizes the dimensionless nor- 

mal force as a function of the softness parameter = hydrodynamic pressure ^ surface deflection 

^ ' elastic stittness gap thickness 

which characterizes the fluid-induced deformation of the interface. The corresponding cases for a 
spherical slider are treated using scaling concepts. 
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I. INTRODUCTION 



The reduction of friction in practical applications has been studied since antiquity. Pic- 
tographs found in Uruk, located in modern day Iraq, have been dated to ca. 3000 B.C. and 
illustrate the transition from sleds to wheels (see Hamrock & Dowson (1981) for a historical 
overview). While this advance certainly reduced friction, further reductions were possible 
upon the introduction of a viscous lubricating fluid in the axle joints. The theoretical under- 
pinings of fluid lubrication in such geometries can be traced back to the work of Reynolds 
(1886), who studied the mechanics of fluid flow through a thin gap using an approximation 
to the Stokes' equations, now known as lubrication theory. Recent efforts in this technolog- 
ically important problem have focused on modiflcations of Reynolds' lubrication theory to 
account for elastohydrodynamic effects (elastic surface deformation due to fluid pressure), 
piezoviscous behavior (lubricant viscosity change due to high pressure), and thermoviscous 
behavior (lubricant viscosity change due to frictional heating) (Dowson & Higginson 1959; 
O'Donoghue, Brighton &Hooke 1967; Conway & Lee 1975; Hamrock & Dowson 1981). 

Inspired by a host of applications in physical chemistry, polymer physics and biolubri- 
cation, in this paper we focus on the elastohydrodynamics of soft interfaces, which deform 
easily thereby precluding piezoviscous and thermoviscous effects. There have been a number 
of works in these areas in the context of specific problems such as cartilage biomechanics 
(Grodzinsky, Lipshitz & Glimcher 1978; Mow, Holmes & Lai 1984; Mow & Guo 2002), 
the motion of red blood cells in capillaries (Lighthill 1968; Fitz-Gerald 1969; Tozeren & 
Skalak 1978; Secomb et al. 1986; Damiano et al. 1996; Secomb, Hsu & Pries 1998; Feng 
& Weinbaum 2000; Weinbaum et al. 2003), the elastohydrodynamics of rubber-hke elas- 
tomers (Tanner 1966; Martin et al. 2002), polymer brushes (Klein, Perahia & Warburg 
1991; Sekimoto & Leibler 1993) and vesicles (Abkarian, Lartigue & Viallat 2002; Beaucourt, 
Biben & Misbah 2004). Another related phenomenon is that of a bubble rising slowly near a 
wall; the bubble's surface deformation then leads to a lift force (Leal 1980; Takemura et al. 
2002). Instead of focusing on specific applications, here we address a slightly different set of 
questions: How can one generate lift between soft sliding surfaces to increase separation and 
reduce wear? What is the role of geometry in determining the behavior of such systems? 
How do material properties infiuence the elastohydrodynamics? Are there optimal combi- 
nations of the geometry and material properties that maximize the lift force? And finally. 



2 



can the study of soft lubrication lead to improved engineering designs and be of relevance 
to real biological systems? 

To address some of these questions we start with the simple case of two fluid-lubricated 
rigid non-conforming surfaces sliding past one another at a velocity V as shown in figure 
n The viscous stresses and pressure gradient due to flow in the narrow contact zone are 
dominant. For a Newtonian fluid, the Stokes equations (valid in the gap) are reversible 
in time, t, so that the transformation t —>■ —t implies the transformations of the velocity 
V — > —V and the normal force L — > —L. In the vicinity of the contact region non- conforming 
surfaces are symmetric which implies that these flows are identical and therefore L = 0. 
Elastohydrodynamics alters this picture qualitatively. In front of the slider the pressure is 
positive and pushes down the substrate, while behind the slider the pressure is negative 
and pulls up the substrate. As the solid deforms, the symmetry of the gap profile is broken 
leading to a normal force which pushes the cylinder away from the substrate. 

This picture applies naturally to soft interfaces which arise either due to the properties 
of the material involved, as in the case of gels, or the underlying geometry, as in the case 
of thin shells. In §1111 we study the normal-tangential coupling of a non-conforming contact 
coated with a thin compressible elastic layer. If the gap profile prior to elastic deformation 
is parabolic in the vicinity of the contact, the contact is non-conforming. However, if a 
parabolic description prior to deformation is insufficient we refer to the contact as conform- 
ing; e.g. the degenerate case considered in §IVI fVl treats normal-tangential coupling of 
non-conforming contacts coated with a thick compressible elastic layer. In §VII we consider 
the normal-tangential coupling of non- conforming contacts coated with an incompressible 
elastic layer. In Will we treat the normal-tangential coupling of non-conforming contacts 
coated with a thin compressible poroelastic layer which describes a biphasic material com- 
posed of an elastic solid matrix and a viscous fluid (Biot 1941). In Willi we study the 
normal-tangential coupling of a non- conforming contact where one solid is rigid and the 
other is a deformable cylindrical shell. In §1X1 we study a conforming contact: a journal 
bearing coated with a thin compressible elastic layer; Finally, in ^JX]we treat the elastohy- 
drodynamics of 3-dimensional flows using scaling analysis. Figure |21 provides an overview of 
the different geometries and elastic materials considered. 

Our detailed study of a variety of seemingly distinct physical systems allows us clearly 
observe their similarities and to outline a robust set of features we expect to see in any 
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soft contact. In all the cases studied the normal force = contact area ■ characteristic hy- 

drodynamic pressure ■L{rj), where Lij]) is the dimensionless lift and the softness parameter 

hydrodynamic pressure rr , x elastic surface displacement rp i , m j rm 

ri = — — r^ff*- flEreometry) = 1 — ^-r- ■ i . lables 111 and 1111 

' elastic stirtness •' ' characteristic gap thickness 

summarize our results for 77 ^ 1. Increasing increases the asymmetry of the gap profile 
which results in a repulsive elastohydrodynamic force, i.e. in the generation of lift forces. 
However, increasing 77 also decreases the magnitude of the pressure distribution. The com- 
petition between symmetry breaking, which dominates for small 77, and decreasing pressure, 
dominant at large r], produces an optimal combination of geometric and material parame- 
ters, rimaxi that maximize the dimensionless lift, L. Whether or not the normal force has a 
maximum depends on the control parameter: the normal force increases monotonically with 
the velocity, but has a maximum as a function of the effective elastic modulus of the sys- 
tem. This suggests that a judicious choice of material may aid in the generation of repulsive 
elastohydrodynamic forces thereby reducing friction and wear. 



II. FLUID LUBRICATION THEORY 



We consider a cylinder of radius R moving at a velocity V and rotating with angular 
frequency lo and immersed completely in a fluid of viscosity yU as shown in figure H The 
surfaces are separated by a distance h[x\ the gap profile, where the x— direction is parallel 
to the solid surface and the z— direction is perpendicular to it. We assume that the velocity 
and pressure field are two-dimensional and in the region of contact we use a parabolic 
approximation, valid for all non-conforming contacts, for the shape of the cylindrical surface 
in the absence of any elastic deformation. Then the total gap between the cylinder and the 
solid is given by 

2 

M^) = /^o(i + ^) + ^(x), (1) 

with H begin the additional elastic deformation and the characteristic gap thickness 
in the absence of solid deformation. The size of the contact zone a/2/io-R, characterizes the 
horizontal size over which the lubrication forces are important, consistent with the parabolic 
approximation in ((T)). If /iq <^ -R, the gap Reynolds number Reg = ^^/p- ~ ^^1% ^ = 
Re, the nominal Reynolds number. Then, if Reg -C 1 we can neglect the inertial terms and 
use the lubrication approximation (Reynolds 1886) to describe the hydrodynamics. For a 
2-dimensional velocity field v = [vxix., 2;), ^^(x, z)) and a pressure field z) the fiuid stress 
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tensor is 

(Tf = /i( Vv + Vv^) - pi. (2) 
Stress balance in the fluid, V ■ cr_^ = 0, yields 

= 

= -d^p + fid^^v^. (3) 

Mass conservation implies 

= d^Vx + d^v^. (4) 
The associated boundary conditions are 



Vx\z=~H = -V, 


^^\z=ho+^ 


—uR, 


Vz\z=-H = VdxH, 




—xu, 






0, 



(5) 

where we have chosen to work in a reference frame translating with the cylinder. We make 
the variables dimensionless with the following deflnitions 

x=^/2hoRx, z = hoz, p = pQp= — p, 



^fx{V-ujR) , 



f = Po C/ = (^'f, h = hoh', H = Hq H\ 




{V-u^R)^^v',, Vx = {V-uRK. (6) 

Here, Ho{^, ...) is the characteristic scale of the deflection, where E' is the effective 
elastic modulus of the medium and k are the length scales of the system. The pressure 



scaling follows from (jH)) and the fact that x ~ y/hoR. Then, after dropping the primes the 
dimensionless versions of equations ((HI)-© are 

= d^p, 
= -dxP + d^^Vx, 
= dxVx + dzVz, 

'^x\z=—'nH Z 7; '^x\z=l+x'^ Z 7) 

1— Uj' 1 — OJ' 

T] dxH —2xuj' 



'^z\z=—riH -. ,5 '^z\z=l+x^ 



I-UJ'' -™ 

p\x^oo = 0, p|x-^-oo = (7) 
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Here lo' = ujR/V characterizes the ratio of rolhng to shding, the softness parameter rj = 
Ho/ho characterizes the scale of the elastic deformation relative to the gap thickness, which 
is related to the compliance of the elastic material. The dimensionless version of the fluid 
stress tensor ^ is 

ed^v^ + e^d^Vz -p + 2e'^d^Vz 
where e = \/ Solving ^ gives the Reynolds equation (Batchelor 1967) 



(8) 



= d,{6h + h^d,p), (9) 



subject to 



p(oo) = p(-cx)) = 0. (10) 
Note that u' is scaled away. Here, h is the gap profile given by (Q) in dimensionless terms 

h{x) = + riH{x). (11) 

To close the system we need to determine r]H{x), the elastic response to the hydrodynamic 
forces. This depends on the detailed geometry and constitutive behavior of the cylindrical 
contact. In the following sections we explore various configurations that allow us to explicitly 
calculate rjH, thus allowing us to calculate the normal force on the cylinder 

pdA, (12) 

contact 

and determine the elastohydrodynamic tangential-normal coupling. We note that when 
rj = 0, the contact is symmetric (lllj) so that the form of Q implies that p{—x) = —p{x) 
and L = 0. 



III. ELASTIC 'LUBRICATION' THEORY: DEFORMATION OF A THIN ELAS- 
TIC LAYER 

In our first case, we consider a thin elastic layer of thickness Hi <^ 1^ coating the cylinder 
or the rigid wall or both, all of which are mathematically equivalent (figure E)). We first 
turn our attention to determining the surface deflection of the layer for an arbitrary applied 
traction. Throughout the analysis we assume that the surface deflection Hq <^ Hi so that a 
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linear elastic theory suffices to describe the material response. The stress tensor, cr^, for a 
linearly elastic isotropic material with Lame coefficients G and A is 

cr, = G(Vu + Vu^) + AV ■ u I, (13) 

where u = {uxiUz) is the displacement field, and I is the identity tensor. Stress balance in 
the solid implies 

V ■ cr, = 0. (14) 

We make the equations dimensionless using 



z = Hiz', X = A/2/io-Ra;', = hou'^, = Hqu'^ cr^ = pocr'^. (15) 
We note that the length scale in the direction is the depth of the layer Hi; the length 



scale in the x— direction is the length scale of the hydrodynamic contact zone, Iq — IHqR; 
the displacements and Uz have been scaled with the characteristic gap thickness /iq; and 
the stress has been scaled using the hydrodynamic pressure scale following We take the 
thickness of the solid layer to be small compared to the length scale of the contact zone with 



Hi 

Then, after dropping primes, the dimensionless 2-dimensional form of the stress tensor (|TH|l 



C = ^j2K^ ^ ^"^^ restrict our attention to compressible elastic materials, where G ~ A. 



IS 




dzUz + CdxUx 2^^^^^' + ^2^*^^^^ 



(16) 



Here 



'^~2G + Xho~'^ 2G + X h'J' ^ ^ 

is the softness parameter, a dimensionless number governing the relative size of the surface 
defiection to the undeformed gap thickness. Stress balance ()14|1 yields 

dzzUx + C(l + ^)dxzUz + C^(2 + ^)dxxUx = 0, 

G + A G 

dzzUz + ^ 2G + X ^''^^'' C^ 2G + x ^"^""^' " ^' 

so that to 0(C) the leading order balance is 



dzzUx = 0, 

dzzUz = 0. (19) 
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The normal unit vector to the soft interface is n = {—dxUz\z=o, 1), which in dimensionless 
form is 

n= {-edxUz\z=oA), (20) 
where e = sj^- The balance of normal traction on the solid-fluid interface yields 

(Tf ■ n|^=o = (^s ■ n|^=o, (21) 

so that 

dzUx\z=o = 0, dzUz\z=o = -VP- (22) 

At the interface between the soft film and the rigid substrate, the no slip condition yields 

Uz{x, ~1) = 0, Ux{x,-l) = 0. (23) 

Solving (fnH) . and gives us the displacement of the solid-fluid interface 

Ux{x,0) = 0, Uz{x,0) = -r]p = -r]H{x). (24) 

This linear relationship between the normal displacement and fluid pressure is known as 
the Winkler or 'mattress' elastic foundation model (Johnson 1985). In light of ()24|1 we may 
write the gap profile (fTTjl as 

h = l + x'^ + r]p. (25) 

Equations Q, pO|) and form a closed system for the elastohydrodynamic response of 
a thin elastic layer coating a rigid cylinder. We note that Lighthill (1968) found a similar 
set of equations while studying the flow of a red blood cell through a capillary. However, 
his model's axisymmetry proscribed the existence of a force normal to the flow. 

When ?7 <^ 1 we can employ a perturbation analysis to find the lift force experienced by 
the cylinder. We use an expansion of the form 

p = p(o) +r7pW + 0(7^2), /i = +r7/i(^) + 0(r/2), (26) 

to find 

r]°: a,.{6/i(°) + = 0, (27) 

T]' : d^Qh^'^ + 3[/i(°f + [h^ydxp^'^} = 0, (28) 
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where 

/i(°) = l+x2, /i«=pW, (29) 
subject to the boundary conditions 

p(o) (oo) = (-00) = p(^) (00) = p(^) (-00) = 0. (30) 
Solving (I27|)-(jnni) yields 

p^frl^+^ll^ + oh^). (31) 



so that 



/ 



pdx = — T]. (32) 
8 



In dimensional form the lift force per unit length is 

3V27r pI HiVR _ 3V2tt fi\V - ujR)^ HiR^/'^ 
~ 8 2G + X ^0 ~ ~ 2G + A 



(33) 



as reported in Skotheim & Mahadevan (2004b). The same scaling was found by Sekimoto & 
Leibler (1993), but with a different prefactor owing to a typographical error. When 77 = 0(1), 
the system 0, (fTUj) and ((23) is solved numerically using a continuation method (Doedel et 
al. 2004) with rj as the continuation parameter. In figureEJwe show the pressure distribution 
p{x), and gap h{x) as a function of rj. For <^ 1, p(— fg) =0. As 77 increases, h increases and 
the asymmetric gap profile begins to resemble that of a tilted slider bearing, which is well 
known to generate lift forces. However, an increase in the gap thickness also decreases the 
peak pressure ~ (see figure 13^). The competition between symmetry breaking, 

dominant for ?7 ^ 1, and decreasing pressure, dominant for ?7 ^ 1, produces a maximum 
scaled lift force aX r] = 2.06. In dimensional terms, this implies that the lift as a function 
of the effective modulus 2G + A will have a maximum, however, the lift as a function of the 
relative motion between the two surfaces V — uR increases monotonically (see figure Eb). 
In fact, ()33|) shows that the dimensional lift increases as {V — uR)"^ for 77 <^ 1. 



IV. DEGENERATE CONTACT 



In this section we consider the case where the parabolic approximation in the vicinity of 
the contact breaks down. Since rotation changes the nature of the contact region for such 
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interfaces, we consider only a purely sliding motion with u = 0. We assume that the gap 
thickness is described by 



™2n 

h = hoil + j^-^) + H{x), (34) 

where n = 2,3, ... characterizes the geometric nature of the contact and the contact length 
is Ic ~ {hQB?^~^Y^'^^\ We note that we always focus on symmetric contacts. We make the 
variables dimensionless using the following scalings 

h = hoh', X = l*x' = hpR^'^^x', 

p = p*p' = ^v?—^ p'. (35) 

2n 

In §1111 we have seen that for a thin compressible soft layer the pressure and and surface 
deflection can be linearly related by so that the scale of the deflection is h^rj = 2g+\ - 
To find the scale of the deflection for the degenerate contact described by we replace 
Po with the appropriate pressure scale p* so that 

''o 

and the size of the deformation relative to the gap size is governed by the dimensionless 
group 

' 2G + \ho 2G + X .3-^ • ^ ^ 

"'0 

Then the dimensionless version of the gap thickness profile (jH^ is 

/;, = l + x2" + r/p. (38) 

As in ^1111 we employ a perturbation expansion for the pressure field in the parameter r) 
(^ 1) to solve 0, (fTn|l and ()38j) and find the pressure field and the lift for small t]. This 
yields p = Po + V Pi ^-^d gives the following dimensionless result 

L„=2 = ^^V, ^n=3 = 0.8859 r], (39) 
784v 2 

where we have not shown the pressure distribution due to its unwieldy size. In dimensional 
terms, the normal force reads as 

35l7r p*2 ^^i* 351^ n^V^HiR^/^ 



->n=2 



78AV2 2G + X ho 78AV2 {2G + X)h'J/^' 



(40) 



p*^ Hil* n^V^HiR^/^ , , 

L„=3 = 0.8859 ^^ , , I = 0.8859^^ 41 

2G + X ho {2G + X)hl^^ 
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The results for n = 2 and n = 3 are shown in figures IHl and [7| One noteworthy feature of a 
degenerate contact is that the torque experienced by the shder arises from the fiuid pressure 
rather that the shear force because the normal to the surface no longer passes through the 
center of the object as it would for a cylinder or sphere. The ratio of the torque due to shear 
to the torque due to the pressure is 

shear torque ^VR/ho ./^Oxi-i x 
(— -) n < 1. (42) 

pressure torque pic H 
Hence, the dominant contribution to the torque is due to the pressure and 



For ?7 <^ 1, p = + Oirj) so that 



pxdA. (43) 



6x ^ 3 vr 

n = 2 : po = 777^ — TT^, r 



7(l + x4)2' 14^2' 



11(1 + X6)2' " 11 ^ 



_ _ 6x _ 37r(2m - 3) esc ^ 

V. SOFT SLIDER 

To contrast our result for a thin layer with that for a soft slider, we consider the case where 
the entire cylindrical slider of length / and radius R is made of a soft material with Lame 
coefficients G and A. Equivalently, we could have a rigid slider moving above a soft semi- 
infinite half space. Since the deformation is no longer locally determined by the pressure we 
use a Green's function approach to determine the response to the hydrodynamic pressure. 
Following Davis, Serayssol & Hinch (1986), we use the Green's function for a point force 



on a half space since the scale of the contact zone, Ic ~ yhoR ^ R, the cyhnder radius. 
Following Landau & Lifshitz (1970) we write the deformation at the surface due to a pressure 
field p(x, y) as 

^ A + 2G' /■ p{x',y')dx'dy' 
""^"^'y^ 4nG{\ + G)J ^(a;'-a;)2 + (y'-y)2' ^ ^' 

with X, y as defined in figure |H1 Neglecting end effects so the pressure is a function of x 

only, we integrate over y' to get 



47rG(A+G) 
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I dy' 

y/(x'-x)2+{y'-y)^ 

{x—x')'^ 



dx' 



^ Too log « Pix')dx', (46) 



To make the equations dimensionless we employ the following scalings 



x = \/2h^x', y = ly', h = hQh', 



n - hoH , p - pqp - — p , (47j 



''0 

so that the dimensionless gap thickness (fTT|) now reads 



where 



/CO Y 
dx'p{x') log[- -r] (48) 
[x - x'Y 



Comparing with our softness parameter for a thin section we see that 
'^thin layer _ /tt G{G + X) Hi Hi 



^^^ TnX Z - ^ « 1' (50) 



^soft shder (2^ + VhoR 

i.e. a thin layer is stiffer than a half space made from the same material by the geometric 

2x 
(l+a;2)2 



factor For small r] we write p = Po + VPi where Pq = iil^2\2 as in (jnH). Substituting 



into (jlHl) yields 



h=l + x' + r^^^ + 0{rf). (51) 
1 + x"^ 



To order 77, Q and (fTUI) yield the equations for the perturbation pressure, pi. 

o r/. 2n'?<^ 247rx (x^ - 1), 
11: = a.[(l + xY9.pi + -^^^^^^], 

pi(-oo) =pi(cx)) = 0, (52) 



which has the solution 

27r(2a;2 - 1) 
[l + x 

Hence the dimensionless lift force is 



L = r] I p^dx = r] I —— — ^c?x = — r/. (54) 



In dimensional terms the lift force is 

^ Stc fi\V -ujR)'^{X + 2G)R'^ 



(55) 



8 G(A + G) hi' 

Comparing this expression with that for the case of a thin elastic layer, equation ()33p. 
we see that confinement acts to reduce deformation and hence reduce the lift in the small 
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deflection, 77 <^ 1, regime. When rj = 0(1) we solve (fTUjl and PHj) for p{x) using an 
iterative procedure. First, we guess an initial gap proflle hold and use a shooting algorithm to 
calculate the pressure distribution. The new pressure distribution is then used in equation 
(I48|) with Y = 1000 (corresponding to a very long cylinder) to calculate a new gap profile, 
knew If jl^ioihoid — hnew^dx < 10~^ the calculatiou is stopped, else we set hgid = hnew and 
iterate. The results are shown in figures 1^ and ITUl and not surprisingly they have the same 
quahtative features discussed previously, i.e. for 77 <^ 1, L ~ 77; L shows a maximum at 
7] = 0.25 and decreases when 77 > 0.25. The reasons for this are the same as before, i.e. the 
competing effects of an increase in the gap thickness and the increased asymmetry of the 
contact zone. 

VI. INCOMPRESSIBLE LAYER 

In contrast to compressible layers, an incompressible layer {e.g. one made of an elastomer) 
can deform only via shear. For thin layers, incompressibility leads to a geometric stiffening 
that qualitatively changes the nature of the elastohydrodynamic problem (Johnson 1985). 
To address this problem in the most general case, we use a Green's function approach. The 
constitutive behavior for an incompressible linearly elastic solid is 



where u = {ux,Uy,Uz) is the displacement and Ps is the pressure in the solid. Mechanical 
equilibrium in the solid implies V ■ cr = 0, i.e. 



a = G'(Vu + Vu^) -psi, 



(56) 















-dyPs + GV^Uy, 



(57) 



Incompressibility of the solid implies 



= V ■ u = dxUx + dyUy + dzUz. 



(58) 
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For the Green's function associated with a point force, (Tzz\z=o = ~f^{x)^iy) where 5 is a 
delta function, the boundary conditions are 

Ux = Uy = Uz = at z = —Hi, 

o-xz = axy = at z = 0, 
cTzz = -fS{x)6{y) at z = 0. (59) 

We solve the boundary value problem (jFfj) - (jKn|) by using a 2-D Fourier transforms defined 
as 

/OO /"OO 
/ Uxikx,ky,z)e-'^''^^+''yyUkxdky, ... (60) 

Then equations (IHZI)-(jSni) in Fourier space are 

= ikxPs + G{-klux - klux + dzzUx), 
= ikyPs + G{-kluy - kyiiy + dzzUy), 
= -dzPs + G{-kluz - k^Uz + dzzUz), 

= —ikxUx — ikyUy + S^-u^, (61) 

subject to the boundary conditions 



(^xz = = dzUx — ikxUz at z = 0, 



(^yz = = dzUy — ikyUz at z = 0, 



^zz = - f = ~Ps + 2GdzUz at z = 0, 

Ux = Uy = Uz = at 2; = — i/;. (62) 



Solving (jF)T|) and (jH^ for ^^(-z = 0) we find 



2-Hi 



uzikx, ky, 0) - /^^^^ ^ ^^^^^^^ ^ + loM2Hi,/l^y) ■ ^''^ 

Since this corresponds to a radially symmetric integral kernel we can take q = ^fc^ + ky 
and use the Hankel transform, Uz = ^ Joir l) lUzdq, to find the surface displacement 
(Gladwell 1980): 

uz{r)\z.. = ^"(^^^IT^WT^^^M^^' = 4^^^^^' ^'^^ 
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where r = y^x'^ + y'^ and Jo{r q) is a Bessel function of the first kind. To see the form of 
the surface displacement, we integrate ()64|) numerically and show the dimensionless function 
W{r) in figure ITTk . The surface displacement H for a pressure distribution p{x) is found 
using the Green's function for a point force so that 



H 



p{x*) 
AtcG 



oo L"'0 



Jo(rg) 



2Hiq - sinh 2Hiq 
I + 2Hfq'^ + cosh2Hiq 



dq 



dy* > dx 



(65) 



where r = \/ {x — x*)^ + (y — y*Y- In terms of dimensionless variables with the following 
definitions 



P = PoP, q 



q' 



X 



V2h^' 
y = ^/2h^y\ y* = y^2h^y' 



a/ 2hoRx' , X* = \/2hoRx' 



2/io-Rr', h = hoh' 
fi{V -uR) R 



AtxG 27iG 
equation ()65p may be rewritten (after dropping primes) as 

2(q — sinh 2(q 



(66) 



H 



p{x* 



Jo(rg) 



1 + 2Cq^ + cosh 2Cq 



dq 



dy* > dx 



(67) 



H, 



is the ratio of the 



where r = ^/{x — x*y + {y — y*Y-, and the dimensionless group C, 
layer thickness to the size of the contact zone. To understand the form of (jU7j) we consider 
the response of a line force p{x*) = 6{x*) and show the results graphically in figure ITTb. 
Moving on to the case at hand, that of a parabolic contact of a cylinder of length 21, we 
write the dimensionless gap thickness as 



h = 1 + x'^ + 1] H, where r] 



1 Po 

2\/2tx G \ 



1 /i(V" - uR) R 

"■o 



2ti 



G 



(68) 



For r/ ^ 1 we can expand the pressure p in powers of t] and carry out a perturbation analysis 
as in ^IIII This yields a linear relation between the dimensionless lift, L, and the scale of 
the deformation: L = Gi{Qr], where Gi{() is shown in figure IT^. The dimensional lift per 
unit length is 



G^iOfi^V-uRY R' 
hi 



(69) 



As C 



27r G 

0, we approach the limit of an infinitely thick layer, in which case there is no 



stiffening due to incompressibility, so that C(C) 



8 ' 



which is the result for an infinitely 



thick layer. To study the effects of confinement on a thick layer, i.e. ( ^ 1, we approximate 
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(EZI) 



as 



H ^ - 



p{x* 



I' 

-I' L^o 



Mrq)[l - (2 + 4(q + C^g^^g-aCg]^^ I ^y*^^ 



where we are now integrating y* over the dimensionless length of the cyhnder /' 



I 



(70) 



»1, 



since the interactions for deep layers are not limited by confinement effects. Evaluating ()70p 
yields 



H 



p{x*] 



-I' 



(^2 _^ 4^2)5/2 



dy* > dx*. 



(71) 



Integrating (f7T|) with respect to y* and keeping the leading order terms in /' yields 



H 



2x* 



[l + x- 



*2^2 



jlog 




+ 2 log 




{x — x*y 





{X-X*y +4:C 



4(r2 + 

8c^[{x - x*y + 



dx* 



(72) 



[{x-x*)^ + ACY 

where we have used the leading order pressure to evaluate p{x*) = ^-^^^*2)2 + r]p^^^ ■ We 



(73) 



integrate ()72|1 with respect to x* to find 

247rx(x2 - 1) A7rx[x^ + 2x\l + 6C + 6C^) + (1 + 2C)^(1 + 8C + 24^)] 

{1+X^y + (1 + 2C)2]3 

As in §IIH equations and (jHn|l yields the system of equations to be solved for the pressure 
perturbation p^^^: 



[6H + 3(1 + x^H + (1 + x2)=^9,pW] = 
.(i)r_oo) =pW(oo) = 0. 



(74) 



Solving yields 

_ 27r(2a;2 - 1) n[12C^ - 20C + 30 - a;2(45 - 60C + 36C^) + 45x^] 



(1 + X2)4 2(4(1 + X2)3 

which we integrate with respect to x to find the dimensionless lift 



rjp^^^dx = rjCi = rj 



1 



30 



0{c 



(75) 



(76) 



On the other hand, as C we approach the limit where the contact length is much 
larger than the layer thickness leading to geometric stiffening. Due to solid incompressibility 
the layer may only deform via shear and the dominant term of the strain is Vu ~ Since 
V ■ u = we see that if ~ ^ ^ ~ Balancing the strain energy / ^(Vu)^ dV ^ 
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G Hilc with the work done by the pressure poHlc yields H ~ ^^jj^- so that 77^^00 ~ 

^1^. For a thin incompressible layer the characteristic deflection is reduced by an amount 
~ (^^ so that lim^^ooCi(C) ~ C'^- Ci is displayed in figure IT^. For intermediate 
values of C, we computed the results numerically. 

The nonlinear problem arising for i] = 0(1) is solved using Q, (fTUI) . (jH7|) and We 
first guess an initial gap profile hold, then a shooting algorithm is employed to calculate 
the pressure distribution. The new pressure distribution is then used in (p7j) and (jUHj) to 
calculate a new gap profile, knew If J^iQihoid — hnew^dx < 10^^ then the calculation is 
stopped, else we iterate with hoid = knew For C = 1 results are shown in figure IT^- 
d. Not surprisingly, we find the same qualitative features discussed previously: L has a 
maximum due to the competing effects of an increase in the gap thickness and the increased 
asymmetry of the contact zone. 



VII. POROELASTIC LAYER 



Motivated in part by applications to the mechanics of cartilagenous joints, we now turn to 
the case of a cylinder moving above a fluid flUed gel layer. This entails a different model for 
the constitutive behavior of the gel accounting for both the deformation of an elastic network 
and the fluid flowing through it. To describe the mechanical properties of a fluid flUed gel 
we use poroelasticity, the continuum description of a material composed of an elastic solid 
skeleton and an interstitial fluid (Biot (1941); for a review of the literature see Cederbaum, 
Li & Schulgasser (2000) or Wang (2000)). Our choice of poroelasticity to model the gel is 
motivated by the following scaling argument (Skotheim & Mahadevan 2004a). Let V and 
Vip denote gradients on the system scale and the pore scale respectively; Pg is the pressure 
varying on the system scale due to boundary conditions driving the flow, while Pp is the 
pressure varying on the microscopic scale due to pore geometry. Fluid stress balance on the 
pore scale implies that the sum of the macroscopic pressure gradient driving the flow, Vpg, 
and the microscopic pressure gradient, Vi^Pp, is balanced by the viscous resistance of the 
fluid having viscosity fi and velocity v, /iV^^^v, so that the momentum balance in the fluid 
yields 

^Vj^^v - Vp, - Wi^pp = 0. (77) 
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When the pore scale, Ip, and system size, Hi, are well separated, i.e. Ip/Hi <^ 1, equation 
(f77j) yields the following scaling relations 

Pa-^^^^-Pp, (78) 

from which we conclude that the dominant contribution to the fluid stress tensor comes from 
the pressure. The simplest stress-strain law for the composite medium, proposed by Biot 
(1941), is found by considering the linear superposition of the dominant components of the 
fluid and solid stress tensor. If strains are small, the elastic behaviour of the solid skeleton 
is well characterized by isotropic Hookean elasticity. For a poroelastic material composed of 
a solid skeleton with Lame coefficients G and A when drained and a fluid volume fraction 
a, the stress tensor cr is given by the constitutive equation 

cr = G(Vu + Vu'^) + A V ■ u I - apgl. (79) 

The equations of equilibrium are 

V ■ cr = 0, (80) 

where we have neglected inertial effects. Mass conservation requires that the rate of dilata- 
tion of a solid skeleton having a bulk modulus (3~^ is balanced by the fluid entering the 
material element: 

= PdtPg + dtW ■ U. (81) 

/i 

j3 ^ \ + 2(7/3 since the Lame coefficients A and G are for the composite material and take 
into account the microstructure, while (3~^ is independent of the microstructure; for cartilage 
(3^^ ~ IGPa, while G ~ A ~ IMPa. Equations (f?^. (jHUj) and (|HT|) subject to appropriate 
boundary conditions describe the evolution of displacements u and fluid pressure Pg in a 
poroelastic medium. 

We now calculate the response of a poroelastic gel to an arbitrary time dependent pressure 
distribution before considering the specific case at hand. To make the equations dimension- 
less we use the following scalings 



X = ^/2h^x', z = Hiz', t = Tt' = — -^^^^t', Uz = hou[, = hu^ 

V — ujK 

Pg=PoP'g = f^{V -ujR)J--^p'g, p = poP, CT = pq(t' = ^{V ~ujR)J-j(t' . (82) 

V '^0 V 
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We take the thickness of the layer to be much smaller that the length scale of the contact 
zone, C = ^2hoR ^ ^^'^ consider a compressible material, G ~ A. Then after dropping 
primes, the stress tensor (f7^ becomes 

Here, 



(83) 



'^~2G + Xho~'^ 2G + X hT ' 
is the dimensionless number governing the relative size of the surface deflection to the un- 
deformed gap thickness, i.e. the material compliance. Stress balance (jHUj) yields 



= d^^u^ + C ^ da,zUz - m — — — d^pgj + C — ^ — <9xxMx, 



and continuity (j8H) yields 



{dzzPg + CdxxPg) = dtPg + dt{dzUz + C^xUx)- (86) 



To leading order (jH3jl and (jHUj) reduce to 

dzzUx = 0, 

- avdzPg = 0, (87) 

and 

^dzzPg = d,pg + -^d^zUz. (88) 

The normal unit vector to the soft interface is n = {—dxUz\z=o , 1), which in dimensionless 
form is 

n = {-edxUz\z=o, 1), (89) 
where e = \J^^- The balance of normal traction on the solid-fluid interface yields 

(Tf ■ n\z=o = (Ts ■ n\z=o, (90) 

so that 

dzUx\z=0 = 

dzUz\z=o - ar]Pg\z=o = -rjp. (91) 
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At the interface between the soft film and the underlying rigid substrate the no slip condition 
yields 

-1) = 0, Ux{x, -1) = 0. (92) 

Solving dHZD, dnH) and yields 

Ux = 0, 

-riP = -aVPg + dzUz- (93) 

To calculate the displacement at the surface Uz{x, 0, t) we need to determine the fluid pressure 
in the gel Pg. Using in (jHH|) yields 

9tPg - idzzPg = Sdtp (94) 

where 

T k^2Bh,{2G + X) 1 .... 

^ = Vp = Hf{V-uR)[mG + X) + a] ^= P[2G + X) + a ^^^^^- ^''^ 

Here t = y VME jg ii^q time scale associated with motion over the contact length, 
Tp ~ -jtq- is the time scale associated with stress relaxation via fluid flow across the thickness 
of the gel, and 5 ~ ^ the inverse of the fluid volume fraction. Equation (jMj) corresponds to 
the short time limit discussed by Barry and Holmes (2001). The boundary conditions for 
fl94|) are determined by the fact that at the solid-gel interface there is no fluid flux and at 
the fluid-gel boundary there is no pressure jump so that 

dzPg{x, -1, t) = 0, pg{x, 0, t) = p, (96) 

where both conditions are a consequence of Darcy's law for flow through a porous medium. 

Although there is flux through the gel- fluid interface, the Reynolds equation Q for the 
fluid pressure will remain valid if the fluid flux through the gel is much less than the flux 
through the thin gap. A fluid of viscosity /i flows with a velocity v = {vx,Vz) through a 
porous medium of isotropic permeability k according to Darcy's law, so that v ~ -Vp. 



Hence, the total flux through a porous medium of thickness Hi is Vx dz, which will 



Hi 



scale as ~ ^lih?" • Comparing this with the flux through the thin gap h^V leads to 



the dimensionless group Qr = 4^. If Qr ^ 1 we can neglect the flow through the porous 







medium. For cartilage. Hi ~ 1mm and /c ~ 10 mm , and flow through the porous medium 



20 



can be neglected if ho ^ lOnm. This implies that the Reynolds lubrication approximation 
embodied in Q remains valid in the gap for situations of biological interest. 

In response to forcing, a poroelastic material can behave in three different ways depending 
on the relative magnitude of the time scale of the motion r = Ic/V and the poroelastic time 
scale Tp. If T ^ Tp the fluid in the gel is always in equilibrium with the surrounding fluid 
and a purely elastic theory for the deformation of the gel suffices; if r ~ Tp the gel will 
behave as a material with a 'memory'; if r <^ the fluid has no time to move relative to the 
matrix and the poroelastic material will again behave as a solid albeit with a higher elastic 
modulus. In the physiological case of a cartilage layer in a rotational joint the poroelastic 
time scale for bovine articular cartilage is reported to be Tp ^ 20 seconds by Grodzinsky, 
Lipshitz & Glimcher (1978), and Tp ^ 500 seconds by Mow, Holmes & Lai (1984). For time 
scales on the order of 1 second, the cartilage should behave as a solid, but with an elastic 
modulus greater then that measured by equilibrium studies. We consider three different 
cases corresponding to: 

A. Low speed: ^ Tp 

When the cylinder moves slowly, ^ Tp, the time scale of the motion is much larger 
than the time scale over which the pressure diffuses {i.e. 7 ^ 1 ~ 6) so that (j94p becomes 

d,,pg = (97) 

Solving (jnZ|) subject to (jHEI) yields 

Pg = P, (98) 

i.e. at low speeds the fluid pressure in the gel is the same as the fluid pressure outside the 
gel. Equations (|^ and (|^ can be solved to yield 

u,{x,0,t) = -{1- a)7]p{x,t). (99) 

We see that this limit gives a local relationship between the displacement of the gel surface 
and the fluid pressure in the gap, exactly as in the case of a purely elastic layer treated in 

m 
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B. High speed: <ti Tp 

When the time scale of the motion is much smaller than the time scale over which the 
pressure diffuses, i.e. 7 ^ 1 ~ 5, becomes 

dtPg = 5dtP. (100) 

Since the gel is at equilibrium with the external fluid before the cylinder passes over it, 
Pg{x, z, —00) = p{x, —00) = 0, and equation (jlOOp yields 

Pgix,z,t) =6p{x,t) (101) 

Inserting (jl01|) into and integrating yields 

In this limit the fluid has no time to flow through the pores and the only compression is 
due to bulk compressibility of the composite gel, which now behaves much more rigidly. 
The effective elastic modulus of the solid layer is now G^ff ~ ~ IGPa rather than 
Geff ~ IMPa. However, the relationship between the pressure and displacement remains 
local as in ^IIII 

We note that if V ■ u = 0, has no forcing term and Pg = 0. Poroelastic theory does 
not take into account shear deformations since these involve no local change in fluid volume 
fraction in the gel. In this case all the load will be borne by the elastic skeleton. However, 
shear deformation in a thin layer will involve geometric stiffening due to incompressibility 
so that the effective modulus will be G^jf ~ ^^G (Skotheim & Mahadevan 2004b). Hence, 
if ^ 1 the deformation should be treated as an incompressible layer as in WII If 

^^GjS ^ 1, the layer should be treated as in §1111 with an effective modulus 

C. Intermediate speeds: ~ Tp 

When 5 ~ 7 ~ 1 rewriting for the difference between the fluid pressure inside and 
outside the gel, Pg{x^ 2;, t) — p{x, t), yields 

9tiPg -p)- ldzz{pg -p) = (S - l)dtp, (103) 
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with the boundary conditions 



dz{Pg-p){x, -l,t) = 0, {pg-p){x,0,t) = 0. 
We expand Pg — p in terms of the solution of the homogeneous part of ()103|) - ()104|) : 

Pg-P = ^ sin7r(n + -)z. 

n=0 

Inserting the expansion into ()103|) we find 

J2[dtAn + J7i\n + h^An] sin7r(n + = {6- l)dtp. 



n=0 



Multiplying ()106|) with sin7r(m + and integrating over the thickness yields 



TzdtP- 



Solving (Hmj) for An{t) yields 

An{t) = 



7r(n+ 2) J-oc 



Substituting ()1U8|) into ()105p yields the fiuid pressure in the gel 
P9=P + Y1 



n=0 



Tx{n + \) 



Finally, (jHl, © and ^Bj yield 



Uz{l) = V 



(104) 



(105) 



(106) 



(107) 



(108) 



(109) 



(110) 



Since the higher order diffusive modes {n > 0) decay more rapidly than the leading order 
diffusive mode (n = 0), a good approximation to ()110|) is 



—rjH^x) = Uz{0) = 1] 



(111) 



This approximation is similar to that used in Skotheim & Mahadevan (2004a). To simplify 
(llllj) for the case of interest we define 

e= r e^^'-''^dt^pdt'. (112) 
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so that 

dti = -^-^i + dtp. (113) 
In the reference frame of the steadily moving cyhnder ^ (x, t) = ^{x — t) so that 

d.i=^-^i + d.P- (114) 

Integrating the above yields 

l-OO _ 2 

i = - e^^'^'-^'^d^^pdx' (115) 



Consequently, the distance between the gel and cylinder is found from (fTTj) . ()111|) . ()112|) and 
(fTTH|l to be 

h{x) = 1 + + r]{l - a)[p + ^ e"^("'-")9,.pdx'] (116) 



7r 

where 



poff; _ V2RHi^i{V - uR) 
^ (2G + A)/^o hf{2G + X) ' ^ ^ 

where we are considering the case where the bulk modulus of the skeletal material is much 

larger than the modulus of the elastic matrix PG <^ 1, so that implies 6 ^ j^. This 

leaves us a system of equations ©, (fTIHl and pi(i|l for the pressure with 2 parameters: 

7] characterizes the deformation (softness); and 7 is the ratio of translational to diffusive 

timescales. The two limits 7 <^ 1 and 7 ^ 1 of ()116p can both be treated using asymptotic 

methods. For 7^1, (jll6p yields 

h = l + x^ + ri{l-a)p, (118) 



and we recover the limit of a thin compressible elastic layer treated in ^Illl with i] (l — a)r]. 
For 7 < 1, dTTI)!) yields 

h = l + x^ + {l--^){l-a)r]p, (119) 



7r2 

which is the result for a thin compressible layer with 77 — > (1 — ^)(1 — a)?]. When ^ 7 <^ 1 
we expand the pressure field as in writing p = Po + rjpi, where Pq = (j^^^i-^i as in 
Inserting this expression into ()116p yields 

/i = l + x2 + r7(l-a) |j9o + ^^ [l + ^{x-x')]d,,podx'^ + 0{-iT]) (120) 

which can be integrated to give 

' +T— — ^ . fl21) 



. ,^2 

/i = 1 + + 77(1 - a 



7r2(l+x2)2 (1 + X2)_ 



24 



We see that increasing 7 increases the gap thickness and lowers the pressure without in- 
creasing the asymmetry, thus decreasing the hft. In the small deflection limit, rj <^ 1, the 
dimensionless lift force is L = Cp{'j)ri and the lift force in dimensional terms is 

, 2t/2 tt p3/2 

i-C,W(l-a.)^^. (122) 

where Cp is a function of 7 and shown graphically in figure IT^. 

When T] = 0(1), we use a numerical method to solve Q, (fTUI) and ()116p on a finite 
domain using the continuation software AUTO (Doedel et al. 2004) with rj and 7 as the 
continuation parameters. The initial solution from which the continuation begins is with 

= 7 = 0, corresponding to /i = 1 + x^, and p = -jj+^ry- The form of the lift force as 
a function of rj, L{ri,'-f) for various 7 can be almost perfectly collapsed onto a single curve 
after appropriately scaling the ri, L axes using the position of the maximum: i.e. 
where Lmax and rimax are shown in figure IT^.c. 



VIII. ELASTIC SHELL 



For elastic layers attached to a rigid substrate the effective stiffness increases with decreas- 
ing thickness. However, for free elastic shells the effective stiffness increases with increasing 
thickness. To see the effects of this type of geometry in biolubrication problems, we turn 
our attention to a configuration in which a surface is rendered soft through its geometry 
rather than its elastic moduli. We consider a half-cylindrical elastic shell of thickness hg and 
radius R moving with constant velocity V parallel to the rigid substrate, while completely 
immersed in fluid of viscosity /i, as shown in figure El The shell is clamped at its edges, 
which are at a height R + Hq above the rigid solid. The shape of the elastic half-cylinder 
is governed by the elastica equation (for the history of the elastica equation as well as its 
derivation see Love 1944) for 9{s), the angle between the horizontal and the tangent vector 
where s is the arc length coordinate. Balancing torques about the point O in figure ITHl 
gives: M{s + ds) — M{s) — ds sin 9 + dsn^ cos 9 + ^{p cos 9 — fxdzU sin 9 — pdgh sin 9) = 
0. In the limit ds — 0, d^M — nj-sin^ + Uz cos 9 = 0. The x— force balance is 
n^i^s + ds) — n^i^s) — dsipdgh + fid^u) = 0, which as ds ^ yields dgn^ = fid^u + pdsh. The 

force balance is nz{s + ds) — nz{s) +pds = 0, which as ds ^ yields dsUz = —p. We note 
that external forces are applied in the contact region where derivatives taken with respect 
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to X are interchangeable with those taken with respect to s, i.e. d^h ~ dsh and dxP ~ dsP, 
which allows for a consistent framework for the fluid and solid equilibrium equations. This 
yields 

G{X + G)hl 



3(A + 2/i) ^ z , 

where the stress resultants determined by the equations 

= dsUx + ex ■ CTj ■ n, = dgUz + • cr ^ ■ n 



I.e. 



n 2 



(123) 



(124) 



(125) 



where, n = {dgh, —1), and ctj is given by (j21). Then, the pressure in the fluid is governed by 
the Reynolds equation (jH)) 



dssP 



-dsh 



(126) 



Finally, since cylindrical deformations are inextensional, we must complement (|123p - (|126p 
with the kinematic equations 



dgX = cos 6, dsh = sin 6. 



(127) 



where the position of the surface of the elastic cylinder is (X, h). The equations ()123|) - ()127|) 
are made dimensionless with the following scalings 



P = PoP = P , s = -kRs , h = RZ, 



h 



3/2 




X = RX' Tlx = pohon'^, Uz = po^/2hoRn'^. (128) 
After dropping primes, the dimensionless forms of equations (|123j) - (jl27p are written as 

dssP = =^[sV2n{^f' + 3Z'dsp], (129) 



dsU^ 

dsX = IT COS 6 
d^Z = TTsin^. 



R 



dssO = 3v27r -j:r{nxsm9 - \ l —riz cost 

G{X + G) hlhT 



I2R 
ho 



(130) 

(131) 
(132) 
(133) 
(134) 
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To find the scale of tlie elastic deformation ^ = wfiere the maximum displacement of 
the cylinder is of order Hq, we note that the change in curvature is of order ^ so that 
the bending strain is e = ^^^^ and the elastic energy per unit length therefore scales as 
/ Ge'^dA ~ / G{^^ fdA ~ The work done by the fluid is due to a localized torque 

and scales as J p x dx ^ Pq hoR, acting through an angle ~ Balancing the work done 
by the fluid torque po hoHo with the elastic energy — ^-2- yields 

l^m. (135) 

SO that (11341) can be written as 



dssO = -^{n^ sin 9 — ■^'^riz cos 9). (136) 

The system dHSl) for 9 has two dimensionless parameters, r] = 

3v^7r^ ^^|^^^^'* ^'^^3/2 and and is subject to 8 boundary conditions 



p(0) = p(l) = 0, 
X(0) = -1, X(l) = 1, 
Z{0) = Z(l) = l + ^, 



^(0) = ^, ^(1) = (137) 

The first two are a consequence of lubrication theory, while the rest are kinematic boundary 
conditions on the cylinder's lateral edges, which are assumed to be clamped. We note that if 
the cylinder has a natural curvature this will not effect the system of equations p23|) - p27p . 
but may change the boundary conditions. 

We solve p29|) - ()133j) . p36|l and p37|l numerically using the continuation software AUTO 
2000 (Doedel et al. 2004) with rj as the continuation parameter. For < 1, the dimensionless 
lift L = Cgf], where the constant Cs{^) is shown in figure IT7b. However, we see that even 
for rj < 100 the linear relationship remains valid. The dimensional lift force for rj < 100 is 

Figure ^1 shows the gap thickness profile and pressure distribution and figure IT7k shows 
the dimensionless lift for ^ = 0.001. As r] increases the elastic deformation breaks the 
symmetry of the gap thickness profile and results in an asymmetric pressure distribution 
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and corresponding lift. However, the concomitant increase in gap thickness decreases the 
magnitude of the pressure. As for the previous systems considered, the competition between 
symmetry breaking (dominant for small rj) and decreasing pressure (dominant for large 77) 
produces a maximum in the lift. The form of the lift force as a function of i], L(r], ^) for 
various ^ can be almost perfectly collapsed onto a single curve after appropriately scaling 
the ri,L axes; i.e. '^^^^^"""'i^^ '^^ where Lmax and rjmax are shown in figure ITTb.d. 

IX. JOURNAL BEARING 

So far, with the exception of ^IVl we have dealt only with non- conforming contacts. In 
this section we consider an elastohydrodynamic journal bearing: a geometry consisting of a 
cylinder rotating within a larger cylinder that is coated with a soft solid. The journal bear- 
ing is a conforming contact and is a better representation of bio-lubrication in mammalian 
joints in which synovial fluid lubricates bone coated with thin soft cartilage layers. Previous 
analyses of the elastohydrodynamic journal bearing have focused on situations where fluid 
cavitation needs to be taken into account (O'Donoghue, Brighton & Hooke 1967; Conway 
& Lee 1975). As before, we restrict our attention the case where the surface deforms appre- 
ciably before the cavitation threshold is reached so that the gap remains fully flooded. A 
schematic diagram is shown in figure ^1 

We take the center of the inner cylinder, Oi, to be the origin; the center of the outer 
cylinder, Oq, is located at x = e^, z = e^. The inner cylinder of radius R rotates with angular 
velocity uj; the stationary outer cylinder of radius R + + Hi is coated with a soft solid of 
thickness Hi <^ R and Lame coefficients A and G. Here, is the average distance between 
the inner cylinder and the soft solid. Following Leal (1992), we use cartesian coordinates to 
describe the eccentric geometry and applied forces, but use polar coordinates to describe the 
fluid motion. When ^ -C 1, the lubrication approximation reduces the Stokes equations in 
a cylindrical geometry for a fluid of viscosity /i, pressure p, and velocity field v = {vr,vg), 
to (Leal 1992) 

drP = 0, 

^dgp = lldrrVg, (139) 
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subject to the boundary conditions 

Vr = 0, Vg = —ujR at r = R, 
Vr = 0, Vg = at r = R + h, 

p{0)=p{2n). (140) 

Since ^ <^ 1, the continuity equation simphfies to (Leal 1992) 

RdrVr + dgVg = 0, (141) 

and the gap thickness profile simplifies to 

h{9) = ho + e^ cos ^ + sin ^ + H{9), (142) 



where H{6) is the elastic interface displacement due to the fluid forces. As in 9IIT 



so that (fra^ yields 



me)^§f^, (143) 



h = ho + e^cose + e, sin + (144) 

2Gr + A 



Using the following primed dimensionless variables, 



Vr = ujhov'r, Vg = ujRv'g, T = h^r' , h = hoh' 



p = p*p' = —^p', = hoe'^, = hoe'^, (145) 



'-0 

we write ()139|) - ()14H) . ()144p . after dropping the primes, as 



subject to 



drp = 0, 

dgp = drrVg, 
drVr + dgVg = 0, 

= 1 + COS 6* + £2 sin 6' + r^p, (146) 



R 

Vr = 0, Vg = —1 at r = — , 

ho 



R , 

Vr = 0, Vg = at r = \- h, 

ho 



p(0)=p(27r), (147) 
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where, 

' ho{2G + X) {2G + X)hl' ^ ' 

is the softness parameter. As in ^jlll we use fll4(ip to derive the system of equations for the 
fluid pressure 

de{&h + eOep) = 0, 
h = 1 + ex cos 6* + e^; sin 6* + ?7 p, 

p{<S) =p{2'k). (149) 

In addition, fluid incompressibility implies that the average deflection must vanish: 

/»27r /'27r 

/ hde = 2n ^ I pde = 0. (150) 
Jo Jo 

The forces on the inner cylinder are 

I psmede = L^, / pcosede = L^. (151) 
Jo Jo 

Here, Lz is the vertical force and is the horizontal force. We begin with the classical 
solution for a rigid journal bearing (Leal 1992): = 1, = r] = ez = 0, ex = 0.053. 
The following brief symmetry argument shows that = when rj = Lx = 0. Since Stokes 
equations for viscous flow are reversible in time, the transformation u —u implies that 
ez —ez- However, due to symmetry we expect the solution to be a reflection about the 
2;- axis and we conclude that = 0. 

As in previous sections we investigate how elastohydrodynamics alters this picture by 
specifying the eccentricity (e^, ez) and calculating the forces {Lx, Lz) as a function of the 
softness parameter 77. Solutions to (jl49j) - ()15Hl are computed numerically using the con- 
tinuation software AUTO 2000 (Doedel et al. 2004) with t] as the continuation parameter 
and the solution for = as the initial guess. Just as for different geometries analyzed in 
previous sections, the deflection of the surface of the soft solid breaks the symmetry and 
leads to the generation of a horizontal force in the x-direction: > 0. For small deforma- 
tions {rj <^ 1) the dimensionless horizontal force L = Cj{ex)f], where the coefficient Cj{ex) 
is shown in figure CHI In dimensional terms, the horizontal force per unit length for small 
deformations is 

= ^^^'^ \2G + X)ho ' ^ = ^^^'^^WTWo ^ ^ 
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For nearly concentric cylinders, -C 1, Cj{ex) = ll5ex + O(e^). For large eccentricities, 
ex 1, the lubrication pressure diverges and Cj(e^) ~ 12.3(1 — e^-)^'^- For r/ ^ 1 we show 
L, h and p in figure QUI 



X. 3-DIMENSIONAL LUBRICATION FLOW 



The analysis of the 3-dimensional problem of a sphere moving close to a soft substrate 
is considerably more involved. Stone et al. {in preparation) are currently engaged in using 
perturbation methods to calculate the elastohydrodjTiamic lift for the case of a sphere trans- 
lating above a thin elastic layer. Here, we restrict ourselves to the use scaling arguments to 
generalize the quantitative results of previous sections to spherical sliders. The results are 
tabulated in Table 2. In the fiuid layer separating the solids, balancing the pressure gradient 
with the viscous stresses yields 

(153) 

where Ic is the size of the contact zone. Substituting h = ho + Hq with Hq <^ ho we find 
that the lubrication pressure is 



p ~ ~ ^(1 + ^) = ^(1 + n) (154) 



The reversibility of Stokes equations and the symmetry of paraboloidal contacts implies that 
the lift force L = when rj = ^ = 0. For rj <^ 1, we expand the pressure as p = p^^^ + rjp^^K 
Since p = p^'^^ will not generate vertical forces, the the lift on a spherical slider, L^, will scale 

as 

Ls=vl P^'^ dA^rij ^dA ~ r/^. (155) 

To compute L, we need a prescription for the softness rj and the contact radius Ic for each 
configuration. 

(A) For a thin compressible elastic layer ( ^IIIIl . we substitute Ic ~ ^/h^R and rj ~ 
^y^Tf ^ ^f/f into (UnSD to find 
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(B) For a thin elastic layer with a degenerate axisymmetric conforming contact f ^IV|) . Ic 
(/ioi?^"~^)^/^" and r] ~ ^f- ^o that Yields 





2t/2 tj d4-^ I 



2G + A - '15^' 

(C) For a soft spherical slider (or thick layer Hi ^ Ic] Wjl. the deflection is given by so 
that 

A + 2G r p{x',y')dx' dy' pl^ , . 

""^""'y^ 47rG(A + G)] ^(x' - xf + {y' -yf^ G' ^ 



so that 

Plc 



the size of the contact zone Ic ~ y/hoR so that (jlSSj) and ()159|) yield 



(159) 



(D) For an incompressible layer f ^Vl|) we have two cases depending on the thickness of the 
substrate relative to the contact zone characterized by the parameter ( = t^- For C ^ 1; 
Ic ~ y/hoR and r] ~ /^(^-^^) ^ gQ ^j^g^^- p55j) yields 

"■0 

For the case C 1 the proximity of the undeformed substrate substantially stiffens the 
layer. In sharp contrast to a compressible layer, a thin incompressible layer will deform via 
shear with an effective shear strain An incompressible solid must satisfy the continuity 
equation V ■ u = 0, which implies that 7^ ~ Consequently, ^ ~ Balancing the 
elastic energy / G{!^ydV ~ G{^^YHill with the work done by the pressure pHqII yields 



Since U ~ v^VR, (USSl), and dHH) yield 

/i^ iff ^ iJ,^V^HfR 



(163) 



'0 







(E) For a thin poroelastic layer i WH^ . k ~ and r/ ~ ^ TcTf ^ '^f//' so that (fTHK|l 

yields 
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where 7 is the ratio of the poroelastic time scale to the time scale of the motion. (F) For 
a spherical shell slider there are two cases: the thickness of the shell, /i^, is smaller 

than the gap thickness, i.e. hg <^ ho and all the elastic energy is stored in stretching; or 
hs ^ ^0 and bending and stretching energies are of the same order of magnitude (Landau & 
Lifshitz 1970). For a localized force the deformation will be restricted to a region of area (P. 
The stretching energy per unit area scales as GhgHQ/R"^, while the bending energy scales as 
Gh^gHl/d^. The total elastic energy, f/, of the deformation is then given by 

U^'^.'^, (165) 

which has a minimum at c? = y/hgR. Comparing d with Ic ~ y/hoR we see that the 



hydrodynamic pressure is localized if hg > ho. For a localized force d = y/hgR while for a 
non- localized force d = R. The elastic energy of a localized deformation, Ui, and a non- 
localized deformation, Un, are given by 

Ui = ^^. Un = GKHl (166) 

The moment exerted by the hydrodynamic pressure on the spherical shell slider is 

M^pll (167) 

which is independent of h^. The work done by the moment ()167|1 . which acts through an 
angle AO ~ Ho/d, is 

MMr^pll^. (168) 

Balancing the work done by the fluid ()168|) with the stored elastic energy ()166|) for both 
nonlocal and local deformations yields 

""-GAA- "'-^^FX' ''''' 
SO that ()155|) and ()169p yield the lift force on the sphere for the two cases 

l^'V' h 

(G) For the ball and socket configuration, roughly the 3-dimensional analog of the journal 
bearing f^IX|). 1^ ^ R and rj ~ (2G+X)h''^ ^^^^ horizontal force is given by pSSj) : 



fM^HiR' 
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XI. DISCUSSION 



The various combinations of geometry and material properties in this paper yield 
some simple results of great generality : the elastohydrodynamic interaction between 
soft surfaces immersed in a viscous fluid leads generically to a coupling between tangen- 
tial and normal forces regardless of specific material properties or geometrical configu- 
rations, i.e. a lift force that arises due to the asymmetric fluid pressure deforming the 

soft solid which breaks the symmetry of the gap profile. For small surface deformations, 

surface displacement i j.r, j- • i ^ c ■ ^■ ■ t 

77 = 1 — ^-p — *- Tu^-T <^ 1, the dimensionless normal force is hnear m ri. Increas- 

' characteristic gap thickness ' ' 

ing 77 [i.e. softening the material) increases the asymmetry but decreases the magnitude of 
the pressure. The competition between symmetry breaking, which dominates for small 77, 
and decreasing pressure, which dominates for large 77, produces a maximum in the lift force 
as a function of 77, the material's softness. 

Additional complications such as nonlinearities and anisotropy in both the fluid and solid, 
streaming potentials and current generated stresses (Frank & Grodzinsky 1987a,b) would 
clearly change some of our conclusions. However, the robust nature of the coupling between 
the tangential and normal forces illustrated in this paper should persist and suggests both 
experiments and design principles for soft lubrication. 
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FIG. 1: A solid cylinder moves through a liquid of viscosity above a thin gel layer of thickness Hi 
that covers a rigid solid substrate. The asymmetric pressure distribution pushes down on the gel 
when the fluid pressure in the gap is positive while pulling up the gel when the pressure is negative. 
The asymmetric traction breaks the symmetry of the gap thickness profile, thus giving rise 

to a repulsive force of hydrodynamic origin. The pressure profile and gap thickness shown here are 
calculated for a thin elastic layer for a dimensionless defiection rj = 10. 
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FIG. 2: Overview of the different geometries and elastic materials considered. G and A are the 
Lame coefficients of the linear elastic material, where G/A = corresponds to an incompressible 
material, Ki is the depth of the elastic layer coating a rigid surface, Ic is the contact length, Tp 
is the time scale over which stress relaxes in a poroelastic medium (a material composed of an 
elastic solid skeleton and an interstitial viscous fluid), and r = Icjy is the time scale of the motion. 
^1111 treats normal-tangential coupling of non-conforming contacts coated with a thin compressible 
elastic layer. ^IVI treats normal-tangential coupling of higher order degenerate contacts coated 
with a thin compressible elastic layer. ^] treats normal-tangential coupling of non-conforming 
contacts coated with a thick compressible elastic layer. ^ VII treats normal-tangential coupling of 
non-conforming contacts coated with an incompressible elastic layer. MVlll treats normal-tangential 
coupling of non-conforming contacts coated with a thin compressible poroelastic layer. Wllll treats 
normal-tangential coupling of non-conforming contacts between a rigid solid and a cylindrical shell. 
^1-X.I treats elastohydrodynamic effects due to coating a journal bearing with a thin compressible 
elastic layer. ^3 treats elastohydrodynamic effects for 3-dimensional flows using scaling analysis. 



40 




FIG. 3: Schematic diagrams of mathematically identical configurations. Configuration (a) is treated 
in the text. 
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FIG. 4: (a) Pressure p{x) as a function of rj. (b) Gap thickness profile, h{x) = 1 + + rjp as a 
function of 77. The initially parabolic gap thickness profile is broken and the maximum value of 
the pressure decreases as 77 increases. 
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FIG. 5: (a) Dimensionless lift per unit length, L, plotted against r], the softness parameter. L has a 
maximum at rj = 2.06 which is the result of a competition between symmetry breaking (dominant 
for r) <^ 1) and decreasing pressure (dominant for 77 S> 1) due to increasing the gap thickness. 
For small rj asymptotic analysis yields L = ^Tj, which matches the numerical solution, (b) The 



dimensional lift force, Lift, is quadratic in the velocity for small velocities while being roughly 
linear for large velocities. Vq - 
and Lq is corresponding lift. 



linear for large velocities. Vq + "^0-^ = " /^V, — is the velocity and rate of rotation at which rj = 1, 
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FIG. 7: Dimensionless lift for /i = 1 + x^" + ?7p where n = 1, 2, 3. The curves are similar, however, 
they can not be rescaled to a universal curve. 
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FIG. 8: Schematic diagram of a gel cylinder moving through a liquid over a rigid solid substrate. 
The asymmetric pressure distribution pushes on the gel when the fluid pressure in the gap is 
positive while pulling on the gel when the pressure is negative. This breaks the symmetry of the 
gap thickness profile, h{x), and gives rise to a repulsive force of elastohydrodynamic origin. The 
dashed line in the lower right hand denotes the undeformed location of the gel cylinder. 
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FIG. 9: Gap thickness h and pressure p as a function of ij for a soft cylindrical gel slider. We note 
that while the pressure distribution is localized to the region near the point of closest contact, the 
change in gap thickness is spread out due to the logarithmic nature of the Green's function of a 
line contact: h=l + x^ + r]J dx'p{x') log[ . ^,^2 ]. 
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FIG. 11: Green's function for a point force (a) and a line load (b) acting on an incompressible 

TT 

layer of dimensionless thickness C = ^ = I- 
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FIG. 12: For 77 ^ 1, L = Ci{Qr] where Cj(C) is shown in (a), (b) shows the dimensionless hft 
per unit length, L, as a function of rj for = \/hQR/Hi = 1. (c) and (d) show pressure, p, and 
gap thickness, h, as a function of 77. As the thickness of the layer decreases the presence of the 
undeformed substrate below is increasingly felt and the layer stiffens. In the linear regime a stiffer 
layer results in a smaller deformation and concomitant decrease in lift. 
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FIG. 13: For rj <^ 1, L = Cp{'y)r] where 7 is the ratio of translational to poroelastic time scales, 
and Cp{'y) is shown in (a), (b) rjmax, the value of 77 at which the lift is maximum, plotted against 

L[ 

rimaxd) ^ ^]^g curvcs Can 



7. (c) The maximum lift Lmax fis a function of 7- After scaling L and 7y, 

be almost perfectly collapsed onto a single curve. 
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FIG. 14: Schematic diagram of a half-cylinder of radius R, thickness hg and Lame coefficients 
and A moving at a velocity V while completely immersed in a fluid of viscosity /x. The edges of the 
half-cylinder are clamped at a distance R + hg from the surface of an undeformed solid. 6 denotes 
the angle between the tangent to the surface and the x— axis. {X{s),Z{s)) are the laboratory frame 
coordinates of the half-cylinder as a function of the arc- length coordinate s. 
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FIG. 15: Schematic of the torque and force balance for a bent cylindrical shell of thickness hg and 
Lame coefficients G and A subject to a traction {—jidzU — pdxh,p) apphed by a viscous fluid, x 
and z arc coordinates in the reference frame of the rigid solid, while s is the arc-length coordinate 
in the shell. M{s) = ^^^^^dgsO is the bending moment. 
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FIG. 16: (a),(b) pressure distribution, p{X), as a function of the softness, rj. As rj increases the 
asymmetry of the pressure distribution increases and the maximum pressure decreases. (c),(d) 
shape of the sheet, where X{s) and Z{s) are the coordinates of the center hne in the laboratory 
frame. We see that the point of nearest contact is pulled back and the symmetry of the profile is 
broken by the forces exerted by the fluid on the cylindrical shell. ^ = 10~^. 
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FIG. 17: (a) dimensionless lift for a cylindrical shell for ^ = 10"'^. For 7] < 100. L = Cs{^)ri, 
where Cs is shown in (b). The form of L(r], ^) can be almost perfectly collapsed onto a single 
curve after appropriately scaling the 77, L axes, i.e. |''"^°^(/'o/g)^ where Vmax and L^^ax are shown 
in (c) and (d) respectively. 
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FIG. 18: Schematic diagram of the modified journal bearing geometry in which the larger cylinder 
of radius R + Hi + ho has been coated by a soft solid of thickness Hi having Lame coefficients G 
and A. The larger cylinder's axis is located a distance in the x— direction and a distance in 
the z— direction from the axis of the inner cylinder of radius R. The average gap thickness is ho. 
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FIG. 19: For small ij the dimensionless horizontal force, L = Cj{€x)rj, where the coefficient Cj is 
plotted above. 
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FIG. 20: (a) Dimensionless horizontal force, L, acting on the inner cyhnder as a function of t], a 
measure of the surface deflection; (b) the corresponding gap thickness profiles; (c) the corresponding 
pressure profiles, = 0, e^; = 0.053. 
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Tables 
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TABLE I: Summary of results for small surface deflections. *Uppcr row corresponds to n = 2, 
while the lower row corresponds to n = 3 and the undeformed dimensionless gap thickness profile 
ish = 1 + 



Geometry 


Material 


Surface displacement 


Lift force/unit length 


T^rnn 1 Enrt^T 
A. lilil idVtJi 


V-.'UlliJJl cool Uic 

QGi'ir' Gnliri 
cido vIkj O'JiivJ. 


AT ^,V HiR^^ 
v^2G+A ,,3/2 


2,%/2-K IJ.^V^ HiR^^ 
4 2G+A ^7/2 

"■0 


Thin layer with 

d ge n r a 1 c nt a c t * 


Compressible 

elastic solid 


iG+A h'' 

2G+A , ii (i 
"() 


35l7r M2y2 ^^^9/4 
784V2 2G+A ,^17/4 


Soft slider 


Elastic solid 


1 ^lV{\+2G) R 
2iT G(A+G) ho 


37r2 ^i^V2(A+2G) ij2 
8 G(A+G) /ij^ 


Thickness ~ y^Rho 


Incompressible 
elastic solid 


1 fJ.V R 
27r G ho 


cm ^"^2 r2 
27r G feg 


Thin layer 


Poroelastic 


V^l-L ^)2G+\ ,,3/2 


"■n 


Cylindrical shell 


Elastic solid 


o. /7T„2A'V(A+2G) R7/2 


f5^AT_2^ .A,nMV'(A+2G) R9/2 


•^V^^ G[X+G) ^3^,1/2 


§. M 


Journal bearing 
thin layer 


Elastic solid 


HiiiRl^oj 
/ig(2G+A) 


p / ^im'^R^ H,R3 
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TABLE II: Summary of results for small surface deflections and spherical sliders. 



Geometry 


Material 


Lift force 


Thin layer 


Compressible 
elastic solid 


,,2t/2 rr d2 

h^{2G+X) 


Thin layer with 
degenerate contact 


Compressible 
elastic solid 




2G+X 


O C±_ 1 • 1 

Sort slider 


Elastic solid 


^2y2 ^5/2 

G .5/2 


Thickness > \/RhQ 


Incompressible 
elastic solid 


G ^5/2 


Thickness <C y/Rha 


Incompressible 
elastic solid 


,,2t/2 U p3/2 
fJj V iliJX ' 

G .5/2 


Thin layer 


Poroelastic 


2G+A /^a 


Cylindrical Shell 


Elastic solid 


^2y2 ^4 
G ^5/2^3/2 


hs > ho 






Cylindrical Shell 


Elastic solid 


^5/2 


hs < /lo 






Journal bearing 
thin layer 


elastic solid 


2G+A /ig 
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